Half-magnetization plateaux in Cr spinels 
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Abstract. Magnetization plateaux, visible as anomalies in magnetic susceptibility at low 
températures, are one of the hallmarks of frustrated magnetism. An extremely robust half- 
magnetization plateau is observed in the spinel oxides CdCr2 04 and IIgCr2 04, where it is 
accompanied by a substantial lattice distortion. We give an overview of the présent state 
experiment for CdCr2 04 and HgCr2 04, and show how such a half-magnetization plateau arises 
quite naturally in a simple model of thèse Systems, once coupling to the lattice is taken into 
account. 



1. Introduction 

Spinels AB2X4, where A and B are métal ions and X={0, S, Se. . . }, are among the most 
ubiquitous of crystal structures. While the undistorted spinel lattice has overall cubic symmetry, 
the network of B-ion sites forms a highly frustrated pyrochlore lattice, built of corner-sharing 
tetrahedra. Simple nearest-neighbour interactions on the pyrochlore lattice cannot select a 
unique magnetic or charge ordered ground state Jj|. It is this fact which makes the properties 
of magnetic spinels — and in particular their high field behaviour — so interesting. 

Chromium spinel oxides ACr204 offer the opportunity to study magnetic frustration in the 
absence of charge and orbital effects. The Cr ion has a strong Hund's rule coupling and lives in 
an octahedral crystal field. Therefore, if A={Zn, Cd, Hg. . . } is a divalent métal ion, the Cr will 
be m a [Ar]3cf high -spin state for which ail of its t2g symmetry orbitals are singly occupied, 
and ail of the Cg symmetry orbitals are empty — see Figure ^ The resulting bulk state is a 
Mott insulator in which ail charge and orbital effects are quenched. The magnetic properties of 
a spinel oxide like CdCr204 therefore dépend by the way in which a set of 5 = 3/2 spins behave 
on a pyrochlore lattice. 

This is a far from simple problem. The minimal magnetic model suggested by the chemical 
structure of ACr204 is 
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Figure 1. Structure of spinel oxide ACr204. O ions -^^ ^^^^^ third-neighbour 
shown in yehow define an octahedral environment for the bonds 
Qp3+ ^Qj-^g^ which have a fuh t2g shell of 3d électrons with 
moment S = 3/2. Detailed control of lattice parameters can 
be achieved by substitution of différent A-site métal ions. 



where Ji is a Heisenberg exchange interaction, h the applied magnetic field and the sum {ij)i 
runs over the nearest-neighbour bonds of a pyrochlore lattice (Figure |2]). In conventional 
magnetic insulators, exchange interactions select a unique ground state for which the classical 
energy per bond is minimized. However, because of the spécial corner-sharing geometry of the 
pyrochlore lattice, the classical energy of any state can be rewritten as 
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where M = (Si -|- S2 + S3 -|- S4)/4 is the magnetization per spin of a given tetrahedron and the 
sum runs over ail tetrahedra. The three components of M = {M^, My, M^) fix only three of the 
four independent classical spin angles per tetrahedron. Thus nearest-neighbour interactions do 
not select a unique ground state, but rather an entire manifold of states for which M = h/8Ji 
in every tetrahedron. Cr spinels can therefore be expected to exhibit very exotic magnetic 
properties. 



At high températures, Cr spinels exhibit very little anisotropy in their magnetism. Their 
susceptibility has the form x = C/{T+9cw), where C ^ 1.9 emu/mol K is in the range expected 
for a S = 3/2 moment with gi ^ 2.0. Values of Ocw vary from 370 K in MgCr204 to 32 K 
in HgCr204 (see Table suggesting dominantly antiferromagnetic (AF) exchange interactions 
between Cr spins. However the Curie-Weiss form of the magnetic susceptibility persists down 
to very low températures, and magnetic order is achieved only for T < Tjy <^ Ocw, where it is 
accompanied by a structural transition. In CdCr204, ZnCr204 and MgCr204, this changes the 
overall crystal structure from cubic to tetragonal and introduces a weak easy-plane anisotropy 
in the magnetism; in IIgCr204 the change is from cubic to orthorhombic (Table 



Table 1. Ordering températures and crystal symmetries of Cr spinels in zéro applied field. 
Values of Tjv and 6cw are for samples prepared by Ueda et al. 0. The nearest neighbour 
exchange Ji is estimated from 6cw ~ zJiS{S + l)/3. 
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The thermodynamic properties of ACr204 at high température and in zéro field are typical of 
a wide range of frustrated magnets, and entirely compatible with the minimal model Equation^ 
The spécial feature of thèse Systems is the structural transition accompanying magnetic order at 
low température. This enables the System to lift the infinité classical ground state degeneracy 
exhibited by Equation Q — a mechanism dubbed order by distortion 5 . 

The question addressed in this paper is what happens to the magnetic properties of Cr spinels 
in applied magnetic field. The scale of exchange interactions in the Cd and Hg compounds implies 
the bulk of their magnetization process occurs over a range of fields now routinely accessible 
in large pulsed magnets. As we shall see, the order by distortion mechanism remains active in 
applied magnetic field, and has some very dramatic conséquences. 

2. Magnetization process and half magnetization plateau 

The magnetization process M{h) of a classical Heisenberg AF on the pyrochlore lattice at zéro 
température is linear up to a saturation field of 8Ji — cf. Equation ((2j. Under the assumption 
that longer range interactions are weak, values of Ji ~ 9 K and Ji « 4 K can be inferred from 
Bcw ~ zJiS{S + l)/3 for CdCr204 and HgCr204, respectively (Table P). Thèse suggest that 
the saturation field is experimentally accessible for both compounds. 

Measurements of M{h) for CdCr204 and ZnCr204 were initially made in pulsed magnetic 
fields of up to 48 T at the ultrahigh magnetic field facility of the Institute for Solid State Physics 
(ISSP) in Tokyo [S]. More recently thèse measurements have been repeated for a range of Cr 
spinels, including HgCr204 |2]. Measurements on CdCr204 have been made at fields of up to 
70 T at the High Magnetic Field Laboratory in Osaka 7 , and preliminary measurements on the 
same System at fields of up to 100 T made at ISSP 8 . 

Results for a variety of Cr spinels taken at ISSP are shown in Figure 01 The magnetization of 
ail compounds is indeed linear over a wide range of (low) magnetic fields. However both CdCr204 
and HgCr204 exhibit a broad magnetization plateau for intermediate values of magnetic field. 
The plateau in CdCr204 has its onset at 28 T jôi, and extends up to approximately 60 T 8, 71. 
It has a magnetization of l.b^s per Cr — exactly half the full moment of an S" = 3/2 Cr^+ ion. 
The plateau in HgCr204 extends from 10 T to 27 T, with a saturation field approaching 50 T. 
Once again the plateau occurs for exactly half the saturation magnetization. It is somewhat 
more rounded than in CdCr204; this can be attributed to disorder in the polycrystalline sample 
used [2]. 

The transition into the plateau state is strongly first-order in both cases, displaying 
a marked hystérésis, and is accompanied by a large change in the unit cell volume of 
\/^V\ /V ^ 0.33 X 10-3 in CdCr204 and \ l^V \ /V ^ 1.6 x 10"^ in HgCr204 El- In CdCr204, 
where good single crystals are available, a colossal négative magnétostriction IS.L/ L = —4 x 10"^ 
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Figure 3. Magnetization processes of a 
variety of single crystal and polycrystalline 
samples of Cr spinels taken at 1.8K, show- 
ing dramatic half-magnetization plateaux in 
CdCr204 and HgCr204. Data taken from |2j. 
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Figure 4. Colossal négative magnétostriction 
accompanying first-order transition into half- 
magnetization plateau state in CdCr204. 
Data taken from (B]. 



can be resolved parallel to field in the [111] direction, accompanied by a magnétostriction of 
similar order in the [110] direction (i.e. perpendicular to the field) — see FigurelU This scale of 
magnétostriction would not be unusual in a System with Jahn-Teller active eg orbitals, but is 
exceptional for an oxide in which orbital degrees of freedom are quenched. Taken together thèse, 
results suggest that the magnetization plateau has its origin in a very strong exchange striction. 
Coupling to the lattice appears to be somewhat stronger in HgCr204 than in CdCr204. This 
is also reflected in a higher ratio of T]\f/6cw for the magnetically ordered phase in zéro field 
(Table lU). 

The transition out of the plateau state at high field is of second order in both compounds. 
Above this second critical field, the magnetization process is again approximately linear, 
suggesting that the new phase arises simply through a canting of the ordered moments of the 
plateau phase. For HgCr204, it is possible to map out the entire magnetization process up to 
saturation at /i ~ 45 T — see Figure 01 A further, first-order magnetic phase transition is 
observed at /i ~ 37 T. While no direct structural data measurements presently exist, this is 
presumably also accompanied by a change in crystal structure or bond length. 

As the température of the sample is raised, the magnetization plateau acquires a finite slope. 
Above a characteristic transition température Te ^ 10 K (CdCr204), Te Ri 7 K (HgCr204) the 
plateau washes out entirely and the magnetization process M{h) remains linear up to fields 
approaching saturation. Thermal fluctuations clearly contribute to the stability of the plateau 
state; its width in h increases with increasing température T, and its transition température 
is somewhat higher than that for the (canted) Néel order at low fields. Quantum and thermal 
fluctuations are known to favour collinear spin configurations |j9|, and in the light of thèse 
results it seems natural to attribute the half-magnetization plateau to a collinear "liitiid" spin 
configuration with three up and one down spins per tetrahedron. It is also interesting to note 
that the critical field for the transition from the canted Néel state into the half-magnetization 
plateau is almost independent of température. 

In summary — the magnetic phase diagram for IIgCr204 deduced from pulsed field 
measurements is shown in Figure It exhibits four phases, which we can identify as : 
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Figure 5. Magnetic phase diagram of 
HgCr204 as a function of température, taken 
from P]. Solid lines dénote first order phase 
transitions; dashed hnes second order. 
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Figure 6. Phase diagram of effective spin 
model Equation Q taken from Jl], as a 
function of magnetic field h and dimensionless 
couphng constant b. Sohd hnes dénote first 
order phase transitions; dashed hnes second 
order. Ah order parameters can be classified 
in terms of the irreps of the tetrahedral 
symmetry group T^. 



• a canted Néel state at low field 

• a cohinear half-magnetization plateau state for intermediate field 

• a second canted state approaching saturation 

• a saturated paramagnetic state for large magnetic field 

Récent elastic neutron scattering data suggest that the half-magnetization plateau state 
in HgCr204 has cubic symmetry and a 16-sublattice magnetic order 10,. The saturated 
paramagnetic state at high field is, presumably, also cubic. We return to thèse points below. 

3. Theory of half magnetization plateau 

Magnetization plateaux are very common in frustrated magnets. They are known to occur, 
for example, in both the Ising and Heisenberg jl2| ll.Sj models on the triangular lattice. 
The conventional explanation for magnetization plateaux is that quantum and/or thermal 
fluctuations act to favour a collinear state in which ail spins point either parallel or antiparallel 
to the applied magnetic fleld. This state may (or may not) exhibit long range magnetic order, 
but since the number of spins in the System is an integer, it must exhibit a magnetization which 
is a rational fraction of the full moment — M = 1/3 in the case of the triangular lattice. This 
leads to plateaux which are weak (being a fluctuation effect), and uncorrelated with lattice 
distortion. 

The essential features of the half-magnetization plateau seen in Cr spinels are its extrême 
robustness and the large accompanying lattice distortion. Any model which seeks to explain it 
should therefore include the coupling between spins and the lattice, as well as the interactions 



between spins. We therefore consider the extended Heisenberg model 



T-C = Jii'i- - OiiPij)SiSj + —pij -h 



ES. 



(3) 



where K is an elastic constant, and the bond variable pij measures the change in séparation of 
two Cr ions relative to the state with the lowest elastic energy |14j . 

The physics of the System is most easily demonstrated if we make two further simplifying 
assumptions — i) we treat the spins as classical variables S, = {sin 9i cos (pi, sinOisin (pi, cos9i) 
and ii) we assume that the spins couple to phonons with crystal momentum q = 0, i.e. that 
there is an overall 4-sublattice order. Under thèse assumptions, we can eliminate the lattice 
variables pij from the problem entirely and work with the pure spin model 



where b = Ja^ /2K is a dimensionless measure of spin-lattice coupling. The biquadratic term 
in Equation (0]) is of a form which favours coplanar or collinear spin configurations, and as 
such is known to mimic quantum and thermal fluctuations j^. A more gênerai choice of spin- 
lattice Hamiltonian can lead to more complex interactions in the effective spin model. However 
thèse do not change the essential physics of the problem, and for our purposes it will generally 
be sufîicient to enforce four-sublattice order by adding a ferromagnetic (FM) third-neighbour 
interaction J3 < 



to the effective spin model Equation 0] (see Figure |2j. The conséquences of relaxing the 
assumptions (i) and (ii) will be discussed below. 

The phase diagram and magnetization process of the effective spin model Equation @ at 
T = are shown in Figure El and Figure [7| The model exhibits four dominant phases : 

• a 2:2 coplanar canted state with tetragonal symmetry for low field 

• a 3:1 collinear (uuud) half-magnetization plateau state with trigonal symmetry for 
intermediate field 

• a 3:1 coplanar canted state with trigonal symmetry for fields approaching saturation 

• a saturated (uuuu) state with cubic symmetry for large magnetic field 

Thèse can be classified in terms of the two-dimensional E, three-dimensional T2 and trivial Ai 
irreducible représentations (irreps) of the tetrahedral symmetry group The predicted 

magnetization M{h) for b ~ 0.2 is in striking correspondence with experiments on Cr spinels. 

Furthermore, we can calculate the changes in bond lengths which occur in each of thèse phases, 
and do indeed find a giant magnétostriction on entering the half-magnetization plateau, as seen 
in experiment. Interestingly, the (négative) sign of the magnétostriction seen in CdCr204 implies 
that the strength of AF exchange interactions Ji increases rather than decreases with bond 
length, i.e. a < in Equation (jS)). This is by no means impossible in a complex magnetic oxide 
with competing exchange paths, and appears to be compatible with the sign of the distortion 
seen in the tetragonal ground state [Tïï] . 

The effective spin model Equation|l]is also amenable to low température expansion and Monte 
Carlo simulation techniques |16| I17j . Thermal fluctuations do not signiflcantly alter the phase 
diagram for parameters relevant to experiment, but impose a finite transition température on 
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Figure 7. Magnetization process of effective 
spin model Equation Q taken from JÂ^, for 
a range of values of the dimensionless spin- 
lattice coupling b. 



0.0 0.1 0.2 0.3 0.4 

T 

Figure 8. Magnetic phase diagram of effec- 
tive spin model Equation @ for J3 = —0.05 
and b = 0.1 at finite température, taken 
from . Both T and h are measured in units 
of Ji. 



each of the phases described above. The predicted magnetic phase diagram is shown in Figure|Sl 
Once again the agreement with experiment (Figure [S)) is, for such a simple theory, striking. 

The four-sublattice theory presented above is almost certainly too simple, in the sensé that 
the conditions in real Cr spinels need not favour four-sublattice states. Indeed, where known, 
the Néel and half-magnetization phases of Cr spinels invariably have a more complex structure. 
However the présent theory is noue the less very successful in explaining their magnetization 
process M (h). This apparent pièce of serendipity is in fact a profound conséquence of the 
extrême frustration of the pyrochlore lattice. The theory presented above accurately describes 
the behaviour of individual tetrahedra in applied magnetic field. Thèse tetrahedra may be 
assembled into arbitrarily complex magnetically ordered states without qualitative change in 
their thermodynamic properties (e.g. magnetization). In fact it is even possible to relax the 
condition of long range magnetic order entirely, and obtain a spin liquid state which exhibits 
much the same magnetization process |16) . 

4. Concluding comments 

The magnetization plateau in Cr spinels is a very récent discovery, and there is every reason 
to believe that much more can be learnt about Cr spinels and other related Systems. The 
mechanism proposed to explain the plateau — a coupling to lattice degrees of freedom — is 
very gênerai and can be expected to function in many other frustrated magnets and magnetic 
molécules. 

Many theoretical avenues also remain open, including the examination of quantum effects 
the ab initio calculation of spin and lattice interactions [10], and the extension of symmetry 
analysis and Monte Carlo simulation to more realistic models j2()j . 

Half-magnetization plateaux are only one of the many interesting new conséquences of 
frustration seen in spinel oxides. However the possibility of conducting further experiments 



and making a quantitative comparison with theory means that they remain one of the most 
rewarding problems to study. 
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